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Abstract. In this paper, we propose a preconditioner based on the shift-splitting method for 
generalized saddle point problems with nonsymmetric positive definite (l,l)-block and symmetric 
positive semidefinite (2, 2)-block. The proposed preconditioner is obtained from an basic iterative 
method which is unconditionally convergent. We also present a relaxed version of the proposed 
method. Some numerical experiments are presented to show the effectiveness of the method. 
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1 Introduction 


We consider the solution of the following large and sparse generalized saddle point problem 


Au 


\ 

-B 



(1) 


where A £ R nxra is nonsymmetric positive definite (x T Ax > 0 for all 0 ^ x £ K"), C £ R mxm is 
symmetric positive semidefinite, the matrix B £ R mxra is of full row rank, x, f £ K", y,g £ R m and 
m < n. It can be verified that the system (1) has a unique solution [1, Lemma 1.1]. Saddle point 
problems of the form (1) arise from finite difference or finite element discretization of the Navier- 
Stokes problem (see [2] and references therein). Several iterative method have been presented 
to solve system (1) or some special cases of it in the literature. The main methods have been 
reviewed in [2]. In [1], Benzi and Golub presented the Hermitian and skew-Hermitian splitting 
(HSS) method to solve (1). Since, in general, the HSS method is too slow to be used to solve [1], 
they used the GMRES method in conjunction with the preconditioner extracted from the HSS 
method to solve (1). Recently, when the matrix A is symmetric positive definite, Salkuyeh et al. 
in [6] have presented a stationary iterative method based on the shift-splitting method to solve (1). 
The proposed method naturally serves a preconditioner for the problem (1). More recently, Cao et 
al. in [3] have considered the same iterative method to solve the system (1) when C = 0. In this 
paper, we consider the problem (1) in its general form and investigate the convergence properties 
of the proposed iterative method and the corresponding preconditioner. 
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2 Main results 


For the sake of the simplicity we use the nations used in [6]. Assuming a,/3 > 0, Salkuyeh et al. 
in [6] considered the splitting A — M a ,p — Af a ,p for the saddle point problem (1) with A being 
symmetric positive definite, where 




al + A B t 
-B pi A C 


and J\f a ,p 


i 


a.1- A -B t 
B pi-C 


( 2 ) 


This splitting gives the following basic iterative method (hereafter is denoted by the MGSS iteration 
scheme) 

M a , 0 u^ = N a , p uW + b (3) 

for solving the linear system (1), where is an initial guess. In continuation, we show that the 
symmetry of the matrix A can be omitted. From Eq. (3), we see that the iteration matrix of the 
proposed method is T a ,p = Hence, the method is convergent if and only if the spectral 

radius of T a ^ is less than 1, i.e., p(Y at p) < 1. 

Lemma 1. ( [6, Lemma 1]) Assume that a and [3 are two positive numbers. If X is an eigenvalue 
of the matrix T at p, then A ^ ±1. 

Lemma 2. Let A £ K™ be a nonsymmetric positive definite matrix. Then, $l(x* Ax) > 0, for any 

o/ie€". 


Proof. Let x = r + is, where r, s £ 
simultaneously. On the other hand, 


Obviously, both of the vectors r and s can not be zero 


x* Ax = ( r T — is T )A(r + is) = r T Ar + s T As + ir T (A — A T )s. 
Hence, 5i(x*Aa;) = r T Ar + s T As > 0. 

Theorem 1. Let A be an eigenvalue of the matrix T and a, [3 > 0. Then |A| < 1. 


□ 


Proof. Let u = ( x;y) be an eigenvector corresponding to the eigenvalue A of r ai( g. Then, we have 
Af a .eu = \M a ,pu which is equivalent to 


(al — A)x — B T y = X(al + A)x + XB T y, 
Bx + (ft I — C)y = —XBx + X(f3I + C)y. 


( 4 ) 

( 5 ) 


According to Theorem 1 in [6] we have x ^ 0. 

Without loss of generality it is assumed that ||x|| 2 = 1. Premultiplying both sides of (4) by x* 
yields 

a — x*Ax — ( Bx)* y = A(a||x ||2 + x* Ax) + A (Bx)*y. (6) 

Since A is positive definite, according to Lemma 2 we have 5ft {x* Ax) > 0. If Bx = 0, then Eq. (6) 
implies 

1^1 \a — x*Ax\ ^/(a — 5ft(x*Ar)) 2 + (3(x* Ax)) 2 ^ ^ 

\a + x*Ax\ ^(pt. + 5ft(x*Ax)) 2 + (3(x*Ax)) 2 < 

We now assume that Bx ^ 0. In this case, from Eq. (5) we obtain 

p{X-l) 


Bx = 


-y + Cy. 


A + 1 


( 7 ) 
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Substituting Eq. (7) in (6) yields 


(1 - A)a - (1 + A )x*Ax = (1 + A) (p^—Ly*y + y*Cy 

Letting p = x*Ax , q = y*y, and r = y*Cy , it follows from the latter equation that 

_ _ . , 1 - A 

am + pqui = p + r, with ui = 


1 + A 


Since ct, p, Si(p) > 0 and q, r > 0, form (8) we see that 


Hence, we have 


|A| 


which completes the proof. 




3f ip) + r 
a + (3q 


|1 — u\ / (1-jRM ) 2 + S(o;) 2 

|1 + u\ y (1 + 5i(w)) 2 + 3(w) 2 


( 8 ) 


□ 


Theorem 1 shows that the MGSS method is convergent and therefore it provides the precon¬ 
ditioner Vmgss = A4 a ,i3 for a Krylov subspace method such as GMRES, or its restarted version 
GMRES(m) for solving the saddle point problem (1). Implementation of the method is as described 
in [ ]. We can also use a relaxed version of the MGSS (say RMGSS) preconditioner 

Vrmgss = ( pi + C ) ' 

for the saddle point problem (1). Similar to Theorem 2 in [6] one may discuss about the eigenvalues 
distribution of the coefficient matrix of the preconditioned system. 


3 Numerical experiments 


We consider the steady-state Navier-Stokes equation 


f -i/Au+(u.V)u +Vp = f, 
\ V.u =0, 


in n = [0,1] x [0,1], 


where v > 0. By the IFISS package [4], this problem is linearized by the Picard iteration and then 
discretized by using the stabilized Q1-P0 finite elements (see [5]). The stabilization parameter is 
set to be 0.1. This yields a generalized saddle point problem of the form (1). The right-hand side 
vectors / and g are taken such that x and y are two vectors of all ones. In Table 1, the generic 
properties of the coefficient matrix have been given. 

We use GMRES(30) in conjunction with the preconditioner Vmgss to solve the saddle point 
problem (1). Numerical results are given in Tabic 1. In this table “Iters” and “CPU” stand, respec¬ 
tively, for the number of iterations and the CPU time (in seconds) for the convergence. To show 
the effectiveness of the methods we also give the results of GMRES(30) without preconditioning. 
We use a null vector as an initial guess and the stopping criterion ||6 — Ax^ H 2 < 10 —9 1|£»||2- In 
the implementation of the preconditioner Vmgss (see Algorithm 1 in [6]), we use the Cholesky 
factorization of pi + C and the GMRES(IO) method to solve the inner systems. It is noted that, 
the inner iteration is terminated when the residual norm is reduced by a factor of 10 2 and the 
maximum of the inner iterations is set to be 40. In the MGSS method the parameters a and P 
are set to be 0.01 and 0.001, respectively. As seen, the proposed preconditioner is very effective in 
reducing the number of iterations and CPU times. 
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Table 1: Numerical results for the test problem with v = 1/50. 


GMRES(30) MGSS 


Grid 

m 

n 

Iters 

CPU(s) 

Iters 

CPU(s) 

8x8 

162 

62 

59 

0.08 

6 

0.90 

16 x 16 

578 

256 

115 

0.28 

8 

2.02 

32 x 32 

2178 

1024 

608 

4.95 

12 

3.58 

64 x 64 

8450 

4096 

3554 

110.1 

28 

21.48 
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